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Abstract 

The A n _i Gaudin model with integerable boundaries specified by non-diagonal K- 
matrices is studied. The commuting families of Gaudin operators are diagonalized 
by the algebraic Bethe ansatz method. The eigenvalues and the corresponding Bethe 
ansatz equations are obtained. 
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1 Introduction 



Gaudin type models constitute a particularly important class of one-dimensional many-body 
systems with long-range interactions. They have found applications in many branches of 
fields ranging from condensed matter physics to high energy physics. For example, Gaudin 
models have been used to establish the integrability of the reduced BCS theory of small 
metallic grains pQ El El |I] and the Seiberg-Witten supersymmetric Yang- Mills theory [Hj. 
They have also provided a powerful tool for constructing the solutions to the Knizhnik- 
Zamolodchikov equation [HI El El El EH! of the Wess-Zumino-Novikov-Witten conformal field 
theory. 

Recently Gaudin models with non-trivial boundaries have attracted much interest (HI 
ITT1 H21 E3 UH EESj • So far, attention has largely been concentrated on Gaudin models with 
boundary conditions specified by diagonal K-matrices. In ^2] > the XXZ Gaudin model with 
boundaries given by the non-diagonal K-matrices in [03 HZj was constructed and solved by 
the algebraic Bethe ansatz method. In this paper we generalize the results in jT2] and solve 
the A n _i Gaudin magnets with open boundary conditions corresponding to the non- diagonal 
K-matrices obtained in [T5] . 

This paper is organized as follows. In section 2, we briefly review the inhomogeneous 
trigonometric vertex model with integrable boundaries, which also services as introducing 
our notation and some basic ingredients . In section 3, we construct the generalized Gaudin 
operators associated with non-diagonal K-matrices. The commutativity of these operators 
follows from applying the standard procedure [EH EDI El UJ] to the inhomogeneous 
trigonometric vertex model with off-diagonal boundaries found in [T5j . thus ensuring the 
integrability of the Gaudin magnets. In section 4, we diagonalize the Gaudin operators 
simultaneously by means of the algebraic Bethe ansatz method. This constitutes the main 
new result of this paper. The diagonalization is achieved by means of the technique of the 
"vertex-face" transformation |21j . Section 5 is for conclusions. In the Apendix, we list the 
explicit matrix expressions of the K-matrices corresponding to the n = 3,4 cases. 
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2 Preliminaries: inhomogeneous A^ij open chain 

Let us fix a positive integer n (n > 2) and a generic complex number 77, and J?(«) G 
End(C n ® C n ) be the R-matrix of the A^ij trigonometric vertex model given by 



= KMEaa ®E aa + J2 \R a ap{n)E w ® Epp + R%(u)Ep a <g> E a pj , (2.1) 

a=l a=£f3 

where Eij is the matrix with elements {Eij) l k = Sj^Su. The coefficient functions are 

sin ( n ) e -iv a > 

K P p{u) = { 1, a = /3, (2.2) 



sin(tt) 
sva(u+rj) 

sin ^) e « q, > q 

sin(«+?7) ' " ' 



R%(u) = { 1, a = /3, (2.3) 



sin(r;) 



L af3 

s'm(u+rj) 

The R-matrix satisfies the quantum Yang-Baxter equation (QYBE) 



e~ iu , a<p. 



R 12 (u x - u 2 )R 13 (u 1 - u 3 )R 23 (u 2 - u 3 ) = ^23(^2 - u 3 )R 13 (ui - u 3 )R 12 {ui - tt 2 ), (2.4) 
and the properties [TH] : 

Unitarity : R\ 2 (u) R 2 i(—u) = id, (2.5) 

+ , . i + , x sin (it) sinfti + nr?) 

Crossing-unitarity : R\\(u)Mz l R% (-u - nri)M 2 = — - - id, 2.6 

sm(u + 77) sm(u + nrj — 77) 

Quasi-classical property : Ri2( u )\ri^o — id- (2.7) 

Here R 2 i(u) = P\ 2 Ri 2 {u)Pi 2 with P\ 2 being the usual permutation operator and U denotes 
the transposition in the i-th space, and 77 is the so-called crossing paramter. The crossing 
matrix M is a diagonal n x n matrix with elements 

M Ql3 = M a 6 a p, M a = e~ 2ia \ a = 1, . . . , n. (2.8) 

Here and below we adopt the standard notation: for any matrix A e End(C n ), Aj is an 
embedding operator in the tensor space C n <S> C n <8> • • •, which acts as A on the j-th space 
and as an identity on the other factor spaces; Rij(u) is an embedding operator of R-matrix 
in the tensor space, which acts as an identity on the factor spaces except for the i-th. and 



j-th ones. The quasi-classical properties ()2.7|) of the R-matrix enables one to introduce the 
corresponding classical r-matrix r(u) 

R(u) = id + T) r{u) + 0(r] 2 ), when 77 — > 0, 

r ( u ) = J2{ r ^ E ^® E ^ + r 5^ E ^® E ^} ■ ( 2 - 9 ) 

Here the coefficient functions in the expression of r(u) are given by 



r2(«) = |;M!(«)} 



(2.10) 



r/=0 

One introduces the "row-to-row" monodromy matrix T(u), which is an n x n matrix with 
elements being operators acting on (C™)®^ 

T{u) = R 01 {u + Zl )R m {u + z 2 ) ■ ■ ■ R 0N {u + z N ). (2.11) 

Here {zi\i — 1, . . . , iV} are arbitrary free complex parameters which are usually called inho- 
mogeneous parameters. With the help of the QYBE ()2.4|) . one can show that T{u) satisfies 
the so-called "RLL" relation 

R l2 (u - vfrMTzly) = T^T^Ruiu - v). (2.12) 

Integrable open chains can be constructed as follows Let us introduce a pair of 

K- matrices K~(u) and K + (u). The former satisfies the reflection equation (RE) 

Ri 2 {ui - u 2 )K{(u 1 )R 21 (u 1 + u 2 )K~{u 2 ) 

= K^{u 2 )R 12 { Ul + u 2 )K^{u x )R 21 {u x - u 2 ), (2.13) 

and the latter satisfies the dual RE [23 

R 12 {u 2 - Ul )K+( Ul ) Mf 1 R-zxi-ux -u 2 - nrj) M x K+(u 2 ) 

= M x K+(u 2 )R 12 (-U! -u 2 - nrj) Mf 1 Kf{ Ul )R 2X {u 2 - u x ). (2.14) 

For models with open boundaries, instead of the standard "row-to-row" monodromy 
matrix T(u) ()2.11|) . one needs the "double-row" monodromy matrix T(u) 

T(u) =T(u)K-(u)T- 1 (-u). (2.15) 
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Using ()2.12|) and ([2.130 . one can prove that T(u) satisfies 



#12(1*1 - u 2 )T 1 {ux)Bax{ux + u 2 )Y 2 (u 2 ) = T 2 (u 2 )R 12 (ux + u 2 )T x {ux)R 2 x{ux - u 2 ). (2.16) 

Then the double-row transfer matrix of the inhomogeneous A^lx trigonometric vertex model 
with open boundary is given by 

r(u) =tr(K + (u)T(u)). (2.17) 

The commutativity of the transfer matrices 

[t(u),t(v)] = 0, (2.18) 

follows as a consequence of (J2.4|) - ()2.6j) and ()2.13j) - ()2.14j) . This ensures the integrability of 
the inhomogeneous A^_ x trigonometric vertex models with open boundaries specified by the 
K-matrices. 

3 A n _i Gaudin model with boundaries 

In this paper, we will consider the non-diagonal solutions K ± (u) obtained in ^H], which are 
respectively given by 

n 

K-{u)l = ^^( W )0f( W ;A)^(-«;A), (3.1) 

3=1 
n 

K + (u) s t = ^%( W )0f(-«;A')0f(«;A'). (3.2) 

3=1 

There are n different solutions, each parametrized by an integer I (1 < I < n) and given by 

J 1, 1 < 3 < I, 

r e - 2i ^, 1 <j < I, 

vs 2J> e ^u+ 2 V ) / + 1 < T < n . 

^ sw(£— u— %r)) ' — J — 

Hereafter we choose some particular I without losing generality. In (j3.1[) and (j3.2[) . 0, 
and <p are intertwiners which are specified in section 4. Besides a discrete parameter I, 
the K-matrix K~(u) (resp. K + (u)) depends on continuous parameters £, {A^} (resp. £, 
{Aj}) and p (whose dependence is through the intertwiner-matrix (j4.2[) below). Here and 



throughout, associated with the boundary parameters {Xj} (resp. {A^}) we introduce a 
vector A = J2 1 j=i^j e j ( res P- A' = X]j=i^j e j)' where {ej\j = l,...,n} is the orthonormal 
basis of the vector space C n such that (e^, e^) = 5jk- 

As can be seen below (e.g. (|4.2j) . (|4.11|) and (|4.12|) ). when m is specialized to A, K~(u) 
does not depend on rj. Moreover, the dependence on A n disappears in the final expression 
of K~(u) although it appears in the expression of <pi(u;X). Without loss of generality in 
this paper we will assume A n = 0. Thus the K-matrix K~(u) depends onn + 1 continuous 
free parameters = 1, ... ,n — 1; p, £} for 1 < / < n — 1, and n parameters = 

1, . . . ,n — 1; p] for Z = n. Some explicit matrix forms of the K-matrices K~(u) for the 
n = 3, 4 cases are given as (|A.lj) - (|A.7|) in the Appendix A. 

Let us emphasize that a further restriction 

N 

A' + 7/5>* = A ' ( 3 - 5 ) 

fc=i 

where {j&lA; = 1, . . . , N} are positive integers such that 2 < < n, is necessary for the 
application of the algebraic Bethe ansatz method in section 4. Hereafter, we shall impose 
this constraint. We further restrict the complex parameters £ and £ to be the same, i.e., 

(3-6) 

so that (|3.7|) below is satisfied. Under these two constraits, the K-matrix K + (u) depends 
on the same set of free parameters as the K-matrix K~(u) (see above). Moreover, the 
K-matrices satisfy the following relation thanks to the restrictions ()3.5|) and ()3.6|) 

\im{K + (u) K~(u)} = lim{K + (u)}K~(u) = id. (3.7) 



Let us now introduce the Gaudin operators {Hj\j = 1,2, ■ ■ ■ , N} associated with the 
lomogeneous 
flUD and 



inhomogeneous A^Li trigonometric vertex model with boundaries specified by the K-matrices 



N ( N 

, -1 



I k+j ) 

where ^(u) = ^{K .{u^^Kj {u) , j = !,-■■, N, with = tr {K+(u)R 0j (2u)P 0j }. 

Here {z.,} are the inhomogeneous parameters of the inhomogeneous A^li trigonometric ver- 
tex model and r(u) is the classical r-matrix given by (|2.9|) . 
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The generalized Gaudin operators (|3.8|) are obtained by expanding the double-row trans- 
fer matrix (|2.17jl at the point u = Zj around r\ = 0: 

r(zj) = rizj^o + rjHj + Oir, 2 ), j = l,.--,N, (3.9) 
Kj = ^<Zi)\n=o- (3-10) 

The relations (|2.7|) and ()3.7|) imply that the first term in the expansion ()3.9|) is equal to the 
identity, 

r( % )| r?=0 = id. (3.11) 

Then the commutativity of the transfer matrices {r(zj)} (|2.18J) implies 

[H j} H k ] = 0, j,k = l,---,N. (3.12) 

Thus the Gaudin system defined by (|3.8|) is integrable. Moreover, the fact that the Gaudin 
operators {Hj} (|3.8|) can be expressed in terms of the transfer matrix of the inhomogeneous 
A^}_ 1 trigonometric vertex model with open boundary enables us to exactly diagonalize 
the operators by the algebraic Bethe ansatz method with the help of the "vertex-face" 
correspondence technique, as will be shown in the next section. The aim of this paper is to 
diagonalize the generalized Gaudin operators Hj, j = 1, ■ ■ ■ , N, (|3.8|) . simultaneously. 

4 Eigenvalues and Bethe ansatz equations 

4.1 Intertwining vectors and face- vertex correspondence 

For a vector m G C n , set 

n 

m j = (m,e j ), |m|=^m fc , j = l,...,n. (4.1) 

k=X 

We introduce [TS] n intertwining vectors (intertwiners) {<pj(u;m)\j = l,...,n}. Each 
(j)j(u;m) is an n-component column vector whose a-th component is {(/>[■ (u;m)}. The n 
intertwiners form an n x n the intertwiner-matrix (in which j and a stand for the column 
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and the row indices respectively), with the non- vanishing matrix elements being 

/ e ifi(m) e iF n (m)+p e 2iu \ 

giFl(m) gi/2(m) 

e iF 2 (m) 

■ . _ e ifj (m) 

giFj (m) • . _ 

p ifn~i(m) 



(4.2) 



y e iF n _i(m) g«/n(m) y 

Here p is a complex constant, and {/j(?7T.)|i = 1, . . . , n} and {Fj(m)|j = 1, . . . , n} are linear 



functions of m: 



i_1 1 
fj( m ) = ^2 m k - rrij - -\m\, j = l,...,n, 



fc=i 



Fj(m) 



F n {m) 



1 

5Jm fc --|rn|, j = 1, . . . ,n - 1, 

k=l 

3, 



— m . 

2 1 1 



(4.3) 

(4.4) 
(4.5) 



We remark that as is clear from (|4.2j) . while (p n (u;m) is a function of u, (pj(u;m) does not 
depend on u for 1 < j < n — 1. 

From the above intertwiner-matrix, one may derive the following face-vertex correspon- 
dence relation [THj : 



#12(1*1 - u 2 )(f>i(u 1 ] m) (g> 0j(m 2 ; - rjei) 

= W^f(ui - u 2 )4>k{ui] m - 77e ; ) <g> 0*(« 2 ; m). 



(4.6) 



Here the non- vanishing elements of {W^- (u)} are 
W#(u) = 1, W$(u) 



sin « 



u 



sin (it + 77) 

sin(u+?j) ' ' 

Eg(?Z) e -iu j < k 
sin(u+??) ' ' 



for j ^ k, 



for j 7^ k. 



(4.7) 
(4.8) 



Associated with {Wjj(u)}, one may introduce "face" type R- matrix W(u) 



W(u)= Y,W%{u)E ki ®E lr 

i,j,h,l 



(4.9) 
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Note that the "face" type R-matrix W(u) does not depend on the face type parameter m, 
in contrast to the Z n elliptic case |2Zj, and thus W(u) and R(u) satisfy the same QYBE, i.e. 
W(u) obeys the usual (vertex type) QYBE rather than the dynamical one. Moreover the 
R-matrix W(u) enjoys the quasi-classical property, 

W(u)\ v ^ = id. (4.10) 

For a generic p G C, the determinant of the intertwiner matrix ()4.2|) is non-vanishing and 
thus the inverse of (|4.2|) exists [18] . This fact allows the introducation of other types of 
intertwiners (p and 4> satisfying the following orthogonality conditions: 

£$ a) (u;m) 4 a) (u;m) = (4.11) 

a 

$ Q) K m + net) 0f \u; m + rye,) = 5 tJ . (4.12) 

a 

From these conditions we derive the "completeness" relations: 

Y,4>P(u;m) 4>f {u;m)=8 aP , (4.13) 

k 

^ 4> k a \u; m + r]e k ) % \u\ m + ne k ) = 8 a p. (4.14) 

k 

Corresponding to the vertex type K- matrices ([3.1]) and (|3.2|) . one has the following face 
type K-matrices /C and JC 



= £ #">(«; A - ^(e, - e,)) A), (4.15) 

a,0 

K{\'\u){ = ^ 4>f\- u i A') K+ ( u )p 4>f\u; X' - rj{ej -£»)). (4.16) 

Straightforward calculations show that the face type K-matrices are diagonal 1 

£(A|«)? = ^ £), ^(AK' = ^i(«), i,j = l,...,n, (4.17) 

where functions {^(-u;^) = ^(w)} and {ki(u)} are respectively given by ()3.3|) and ([3.4]) . 

A remark is in order. Although the K-matrices K^lu) given by ([3.1]) and ([3.2]) are 
generally non-diagonal (in the vertex picture), after the face- vertex transformations ([4.15]) 

1 The spectral parameter u and the boundary parameter £ of the reduced double-row monodromy matrices 
constructed from K,{\\u) will be shifted in each step of the nested Bethe ansatz procedure Therefore, 
it is convenient to specify the dependence of K{\\u) on the boundary parameter £ through kj(u; £). 
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and f)4.16|) . the face type counterparts JC(\\u) and JC(\'\u) become diagonal simultaneously . 
This fact enables ones to apply the generalized algebraic Bethe ansatz method to diagonalize 
the transfer matrix r(u) [HOl^JI^] (for n = 2 case, the corresponding transfer matrix was 
diagonalized alternatively by the fusion hierarchy of the transfer matrices with the anisotropy 
value being the roots of unity [3TJ 32J). 

4.2 Algebraic Bethe ansatz 

By means of ()4.13|) . ()4.14p . ()4.16|) and (|4.17jl . the transfer matrix t{u) ()2.17|) can be recasted 
into the following face type form: 

t(u) = K>{*\< T(AuY» = E W ( 4 -!8) 
Here we have introduced the face type double-row monodromy matrix T(m\u), 
T (^'\ u )^ = T{m\u)^\ m=y = 4> v (u;m-r](e^ - e„)) T(u)0 M (-u;m) 

= ^^(M'-^-^TMf^hM')- (4-19) 

a,/3 

Moreover, (|2.1fi|) . (|4.fi|) and (|4.14jl imply the following exchange relations among T(m\u)^ l : 

E E W nl>i - «*») T ( m + < e n + <fc)l«i)£ 
x Wj? a («x + ifc) T(m + r)(e h + e h )\u 2 )% 

= E E T ( m + ^ + c *o)i u »)$ w?*^ + W2 ) 

x T(m + 77(e i2 + e i2 )K)£ I i ! » i - u 2 ). (4.20) 

In the following we will use the standard notation, 

A(m\u) =T(m\u){, Bj(m\u) =T(m\u) 1 j, Cj(m\u) =T{m\u)\_, j = 2,...,n, (4.21) 
T>{{m\u) =T(m\u)i - 5jW(-(2u) A{m\u), i,j = 2,...,n. (4.22) 

In order to apply the algebraic Bethe ansatz method, one needs to construct a pseudo- vacuum 
state (also called a reference state) which is the common eigenstate of the operators A, T> ] - 
and is annihilated by the operators Cj. In contrast to the models with diagonal K ± (u), the 
usual highest-weight state 



1 1 ^ 




fi\ 





® • •• ® 





V : J 
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is no longer the pseudo-vacuum state. However, after the face- vertex transformations ()4.15|) 
and (|4.16p . the face type K-matrices JC(X\u) and JC(X\u) simultaneously become diagonal. 
This suggests that one can translate the trigonometric model with non-diagonal K- 

matrices into the corresponding SOS model with diagonal K-matrices JC(X\u) and JC(X\u) 
given by (|4TH|l - (j4~TH) l. 
Consider the state 



\Q) = fa{-z x - X - (N - 1)7761) ® M-Z2; X — (N — 2)776!) • • • ® 4>i(-z N ; A), (4.23) 

which depends on the boundary parameters {Xj}, but not on the boundary parameter £. It 
is also independent of {zi} because as mentioned before the vector <pi(u; m) does not depend 
on u regardless of m. The state (j4.23|) is the pseudo- vacuum state in the "face picture" since 
it can be shown to satisfy the following equations: 



A(X - Nrje^u) \Q) = h(u; (4.24) 

a sin(2^ 77 
sin(2w + 77) n 2' s 2' 



JV 



tt sin(w + z k ) sin(w - z^) I 
sin^ + ^ + 77) sin(?j — z k +rj) \ 

a,j=2,...,n, (4.25) 
Cj(A-JV?7ei|u)|n) = 0, j = 2, (4.26) 

/3,(A- A^lu) |ft> ^ 0, j = 2,...,n, (4.27) 

as required. 

For later convenience, we introduce a set of non-negative integers {Nj\j = 1, . . . , n — 1} 
with Ni = N and complex parameters {i>jf^| = 1,2,..., Nj + i, j — 0, 1, . . . , n — 2}. As 
in [33J EH E3 121]; the parameters {v^ } will be used to specify the eigenvectors of the 
corresponding reduced transfer matrices (see below). We will also adopt the convention, 

v k = vf\ k = l,2,...,N. (4.28) 



We now apply the generalized algebraic Bethe ansatz method developed in 29J to diag- 
onalize the transfer matrix (j2.17|) with the K-matrices K^^u) given by (j3.1|) - ()3.5|) . We seek 
the common eigenvectors of the transfer matrix of the form 

n 

\v u ...,v N ) = Yl F h ' l2 '-' lN BhW + -r7e 1 |7Ji)/3 i2 (A , + 77(e n + e i2 ) - 2^2) ■ ■ ■ 
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JV-1 



k=l 
N 

'B iN (^' + vJ2^ -Nvei\v N )\Q). (4.29) 



fc=i 



Here the summing indices should obey the constraint, A' + «2~^fcli Q fc = A, where A' and A 
are the boundary parameters which satisfy the restriction ()3.5|) . Taking this constraint into 
account, ()4.29j) may be written as 



\v 1 ,...,v N )= Yl F h ^ N B h (A' + V e h - K (A' + r](e h + e h ) - 2rje 1 \v 2 ) 

i\,...,i N =2 

N-l 

r/ \ I 

"JV — 1 \ » 

k=l 

xB^X-NrjexMlty. (4.30) 



With the help of ()4.18|) . (j4.21|) and ()4.22|) one may rewrite the transfer matrix (j2.17|) in 
terms of the operators A and V 3 - 

n 

t(u) = a {1) (u) A(X'\u) + ~ k f\ u + |) v ( x '\ u ) J r ( 4 - 31 ) 

i=2 

Here we have introduced the function a^(u), 

n 

aW(u) = Y,h(u)W{}(2u), (4.32) 

and the reduced K-matrix ^^^(A'|m) with the elements given by 

iC^(X'\u)i = Sfkf\u), z,j = 2,...,n, (4.33) 

kf\u) = fy(u-|), J = 2,...,n. (4.34) 

Following CHI EI] , we introduce a set of reduced K-matrices {K,^(X'\u)\b — 0, . . . ,n — 1} 
which include the original one )C(X'\u) = JC^(X'\u) and the ones in (j4.33|) and ()4.34|) : 

K {b) (X'\u)i = Sjkf\u), i,j = b+l,...,n, 6 = 0,..., n-l, (4.35) 

~kf\u) = kJu-b^-), j = b+l,...,n, 6 = 0,..., n-l. (4.36) 
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Moreover we introduce a set of functions {a^ b \u)\b = 1, ... ,n — 1} (including the one in 
()4.32j0 related to the reduced K-matrices K.^(X'\u) 

n 

a ®(u) = Wl-(2u) ~kf- 1] (u), 6 = 1, ... ,n. (4.37) 

j=b 

Carrying out the nested Bethe ansatz [2IJ, one finds that, with the coefficients F 4l,i2 '"' ,!jV 
in ()4.3()j) properly chosen, the Bethe state \v\, . . . , vn) is the eigenstate of the transfer matrix 

(Em 

t(u) \v x , . . . , v N ) = A(u; £, {v k }) ■ ■ ■ , vjv), (4.38) 
with the eigenvalue given by 
A(«; £,{«*}) 



^ x sin « + v k + T)) sin w - 7> fc ) 



sin(27i)e" ? f -py sin(w — u fc + 77) sin(w + i> fc + 277) 
sin(2?j + 77) lr~ sin(w — 7^) sin(7t + u fe + 77) 

-pr sin(u + sin(w - z k ) 
svn(u + z k + 77) sin(u - z fc + 77) 



x 



A (1) («+|;e-|, ■ (4.39) 

The eigenvalues {A&'>(i*; £, {7'? ) })|j = 0, . . . , n - 1} (with A(«; f , {u fc }) = A< >(«; f , {^ 0) })) 
of the reduced transfer matrices are given by the following recurrence relations: 

h=1 sm{u + v k ' + rj) sm{u — v k 'j 

I ^3 + 1 / ( 7 ") \ / f ■? 1 \ 

sin(2u)e i?? J j-j- sin(7t — + 77) sin(w + v% + 277) 



sin(2n + 77) J gin ( w _ V C7)) sin ( w + V U) + ^ 

-j-j- sin(7t + z^) sin(w — z^) 



^ sin(?/ + + 77) sin(7i — z% + rj) 

xA0^)( U+ |;^)-f,{4 i+1) }) 
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i = l,...,n-2, (4.40) 

A(n-D( M;e (n-D) = fc(«- 1 )( u )fc n ( W ;^"- 1 )). (4.41) 

The reduced boundary parameters and inhomogeneous parameters } are given by 

= - f, 4' +1) =4 i} + |, j = 0,...,n-2. (4.42) 

Here we have adopted the convention: £ = ^°' ) , z k °^ = z k - The complex parameters {ujjf } 
satisfy the following Bethe ansatz equations: 

sin(2v s + i])e~ tv 



a W (v s )h(v s ;0 : 



sin(2f s + 2rf) 
^ -q sin(w s + ujfc) sin(w s - v k - rj) 



n 



£ =i sin(t; s + w fe + 2?7) sin(t; s - v fc + 77) 
shifty + sin(w s - z k ) 



k^s 

N 



L sin(v s + z k + r]) sin(v s - z k + 77) 



xAW^+Je-?,^}). ( 4 - 43 ) 



a 



2'" 2 
\ « sin(2us + ?])e 



-IT] 



X 



sin(2f?' ) + 277) 

N i+i ■ 1 (i) 1 ■ 1 (i) (;■ 

- - sm(t> s + ) sm{Vs — v k — 77) 

k ^ =1 sm(v ( s j) + v { k j) + 27?) sin(^ i) - v{> + 77) 



n 



k=l sm(vs + 2;^ + 77) sm(vs - + 77) 

xA o- + D ( ,(i) + | e o)_| !{ 4^) })) 

j = l,...,n-2. (4.44) 

4.3 Eigenstates and the corresponding eigenvalues 

The relation (j3.10j) between {Hj} and {t(zj)} and the fact that the first term on the r.h.s. of 
()3.9|) is a c-number enable us to extract the eigenstates of the generalized Gaudin operators 
{Hj} and the corresponding eigenvalues from the results obtained in the previous subsection. 
Introduce the functions 77)} 

^ +1 \u, 77) = 0^ x \u) = a^ +1 \u)k J+1 (u; £- J - V ), j = 0, • • - ,n - 2. (4.45) 
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Then by fl4.40|) - ()4.42|) . the Bethe ansatz equations (j4.43J) and ()4.44|) become, respectively, 

(i) \ _' Nj+i , (j) (j). , (j) (j) v 

sin(2^ ) + 2rj) k ^ =1 sin(^ j) + + 2t?) sm(v { s 3) - v ( k j) + rj) 

= /jc^fow + ri )T\ MvP + vtv + VMvP-vt 1 '-!) 
p u + 2 ; ll sin( ^ + u a-i) + | )sin( ^ 

X jfj sin(^' ) + + f ) sin(t4 j) - v ( k j+1) + 2) ' 

i = 0,---,n-3, (4.46) 

^ (n _i ) ^ i ( w - 2) ^ sin(2t;i"~ 2) + ^e"^ sin(^ n ~ 2) + ^"~ 2) ) sin(^ n ~ 2) - v { k ~ 2) - 77) 

sin(2fi n_2) + 2rj) k ^ S)k=1 sin(v!; n ~ 2) + fj: n ~ 2) + 2rj) sin(i;i n ~ 2) - ^ n ~ 2) + rj) 

= kt 1 \vt- 2) + \)Uv^ + \,(, { "- 1) ) 

Here we have used the convention: v\, — Zk, k = 1, • • • , N. The quasi-classical property 
fETTUJ) of the R-matrix W(u), (IQ7I) and (OHl) lead to the following relations 

pV +1 \u,0) = l, jV" +1 >M) = 0, j = 0,--.,n-2. (4.48) 



Then, one may introduce functions {7^ +1 ^(m)} associated with 77)} 

j = 0,---,n-2. (4.49) 



or] 



r)=0 

Using (|4.48p . we can extract the eigenvalues hj (resp. the corresponding Bethe ansatz equa- 
tions) of the Gaudin operators Hj ()3.8|) from the expansion around 77 = for the first order of 
77 of the eigenvalues ()4.39|) of the transfer matrix r(u = Zj) (resp. the Bethe ansatz equations 
()4.46j) and ()4.47j) ). Finally, the eigenvalues of the generalized Gaudin operators are 

m 

hj = j^(Zj) — {cot(zj + x k ) + cot(zj — x k )} . (4.50) 

k=l 

The parameters {x k ^\ k — 1,2, • • • , Nj+x, j = 0, 1, • • • ,n — 2} (including x k as = xf\ k = 
1, • • • , N\) are determined by the following Bethe ansatz equations: 

Nj+i 

7 (i+i)( x 0')) _ cot (2xi j) ) -i-2 {cot(x^ + x[ j) ) + cot(x^ - xf)} 

k^s,k=l 
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7 0'+2 )(x (i) ) _ {cot(a#> + x ( i~ l) ) + cot(s?'> - x { t 1] )} 



fc=l 

N j+2 



- £ {™t(xi j) + x^ +1) ) + cot(x® - 4' +1) )} > 

i = 0, • ■ ■ , n — 3, (4.51) 
- cot(2^- 2 )) - i - 2 ^ | cot (x^- 2 )+x^ 2) ) + cot(a;(' 1 - 2 )-4"- 2) )} 



■V, 
k=l 



{n - 2) ) - E {^r 2) + 4 n_3) ) + cot(xi- 2 ) - 4- 3) )| . (4.52) 



Here we have used the convention: x k — Zk, k = 1, • • • , iV in ()4.51|) . and the function g( 
is given by 



5 Conclusions 



(4.53) 

r?=0 



We have studied the A n _i Gaudin model with boundaries specified by the non-diagonal 
K- matrices K ± (u), (|3.1|) and ()3.2|) . In addition to the inhomogeneous parameters {zj}, the 
Gaudin operators {Hj}, ()3.8j) . have one discrete parameter 1(1 < I < n) and n + 1 continuous 
free parameters {Xj \j = 1, . . . , n — 1; p, £} for 1 < / < n — 1 (or n continuous free parameters 
{Xj\j — 1, . . . , n — 1; p} for Z = n). As seen from section 4, although the "vertex type" K- 
matrices K^iu) (|3.1j) and (|3.2jl are non-diagonal, the compositions, ()4.15j) and (|4.16)1 . lead to 
the diagonal "face-type" K-matrices after the face-vertex transformation. This has enabled 
us to successfully construct the corresponding pseudo- vacuum state \Vt) (|4.23|) . diagonalize 
the generalized Gaudin operators {Hj} by means of the algebraic Bethe ansatz method, and 
derive the eigenvalues (|4.5Uj) as well as the Bethe ansatz equations (|4.51|) and (|4.52|) of the 
boundary Gaudin model. 
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Appendix: Explicit matrix forms of K-matrices 

In this appendix, we give the explicit matrix expressions of the K-matrices for the cases 
n = 3,4. 

The case: 

There are three types of K-matrices K~(u), which are labelled by the discrete parameters I, 
for the trigonometric A% model. 

• For the case of / = 1, the 7 non- vanishing matrix elements K~{un are given by: 



e 2iu + e p y gin ( w + 



e 2iu + eP y sin ( w + £)J > 



e -2i(X 1+ X 2)+P , 2 . uSin ( u _^) 



K-(u)\ = - s- . : 1 + e 



K-{u)i 
K~(u 



' ~ e 2iu + eP y~ sin(w + f ) J ' 



(e 2iu + eP) sin(u + £) 
2 1 (n sin(u-£) 



e 2iu + e p y s j n ( u _|_ £) 



r u2_ e-^-(»-^sin2^ rW 3_ e - 2 ,^H fAl , 

For the case of / = 2, the 7 non- vanishing matrix elements are given by: 

K~(u){, K~{u)\, K~(u)l are the same as those in (jA. 1|) . 

e 2i(A 1+ A 2 )+ t (n+5) sin2M 

K (u)* = 1, K {u)\ = 



K-(u) 3 2 



(e 2iu + eP) sin(w + f ) 
e 2i\ 2 +i(u+o sin2u 



(e 2iu + eP) sin(w + £) ' 



K-{u)l = ^—(ef-^$.). (A.2) 

v i a e 2iu + eP y S m(w + 0, 
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For the case of I = 3, the 5 non- vanishing matrix elements K (u)^ are given by: 



K-{u)\ 



g2iu _|_ g4iu+p 



2iu 



+ eP 



K-(u)l 



-2iXi+p ^g4iM _ j\ 



3 2m 



+ eP 



K~(u 



p -2i(X l +X 2 )+p / 4m _ i \ 

^, K-(u)l = K-(u)l 



2tu 



1. 



(A.3) 



The case: 



There are four types of K- matrices K (u), which are labelled by the discrete parameters I, 
for the trigonometric model. 

• For the case of / = 1, the 10 non- vanishing matrix elements R'(u)^ are given by: 

K~ 



„2iu 



R~ 
R- 
R- 
R- 
R- 
R- 
R- 



ti = 



U)n = 



U 1 Q = 



g2m gp 

g-2iAi+p 

g2m gp 

; -2i(Ai+A 2 )+p 
„2iii 



sin(w + £) / 5 



1 + e 



1 + e 



,2iu sin ( M -0 

sin(-u + £) J ' 

2m sin (^-0 



M) 4 = 



till = 



«)2 = 



«)2 = 



g^m gp 

e -2j(Ai+A 2 +A 3 )+p 

g2m gp 

e 2i\ 1+ i(u+0 gin 2u 



1 + e 



sin(ti + £)/ 



sin(ti + £) 



^2m _ gp) g in ( u + 



I ( . sin(n — £) 



e 2m — e p ^ ' sin(u + £) 

e -2jA 2 -i(u-0+P s i n 2m 
( e 2iu _ e p) g in ( M + ' 

g-2i(A 2 +A 3 )-j(M-5)+p g j n 2 M 

(e 2iu - e^) sin(u + f) ' 
n4 ,-2» sin (?-«) 



sin(£ + ti) 



(A.4) 



For the case of / = 2, the 10 non-vanishing matrix elements R (u)j are given by: 
R~(u){, R~(u) 2 , R~(u) 3 , R~(u)\ are the same as those in (|A.4j) 



K'{u)l = 1, 



g2i(A 1+ A 2 )+i(«+0 gin 2u 



R-(u)l 



(e 2iu - eP) sin(w + £) ' 
e 2i(\ 2 )+i(u+o sin2ti 



(e 2i » - eP) sin(u + f ) ' 
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sin(u — £) 



3=^ — I" - - 1 " 



e*'" — e p \ sin (it + £) / ' 
e -2iA 8 -i(«-{)+P s i n 2 u ^_ 4 _ _ 2iu sin(£-u) 
(e 2m - e^j sin(w + |) ' ~~ 6 sin(£ + u) 



For the case of / = 3, the 10 non-vanishing matrix elements K (u)j are given by: 

K~(u) 1 , K~(u) 2 , K~(u) 3 , K~(u)\ are the same as those in (|A.4|) 
K-{u)l = K{u)l = 1, 



e 2i(Ai+A 2 +A3)+i(u+g) g j n 2 M 

^ (e 2 ™ - e") sin(u + f) ' 

e 2i(A 2 +A3)+i(«+«) sin2u 

~ ~ (e 2 » - e") sin(u + |) ' 

_ e 2iA 3 +i(«+?) gin 2n 

ii (IT 



3 ( e 2m _ g p) sin ^ + ' 

= __J_(e> + $^|V (A.6) 
• For the case of / = 4, the 7 non- vanishing matrix elements K~(u)j are given by: 

/r»J = —- - — , k-{u)\ = - — -X tl, 

„-2i(Ai+A 2 )+p ( p Uu _ 1 \ 

= " f e^X - ■ 

e -2j(Ai+A 2 +A 3 )+p / g 4w _ -Q 
^( W )4 = e 2i U _ gp > 

K-(^ = K-(^ = K-K = 1. (A.7) 

The above explicit results confirm the following properties for the non-diagonal K-matrices 
of the vertex model: for I = 1, . . . , n — 1, K~(u) depends on n + 1 continuous free 

parameters and has 3n — 2 non- vanishing matrix elements; for I = n, K~(u) depends on 
n continuous free parameters and has 2n — 1 non-vanishing matrix elements. So our K- 
matrices contain more boundary parameters and more non-vanishing matrix elements than 
those found in j^l E3 EH- 
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